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Abstract 

The electroweak phase transition, if proceeding through nucleation and 
growth of bubbles, should generate large scale turbulent flow, which in turn 
generates magnetic turbulence and hence magnetic fields on the scale of tur¬ 
bulent flow. We discuss the seeding of this turbulent field by the motion of 
the dipole charge layers in the phase transition bubble walls, and estimate the 
strength of the produced fields. 

I. INTRODUCTION 

The problem of the generation of magnetic helds is an old one in cosmology The 

general approach is to identify mechanisms for the generation of seed helds B- 0. which 
can become amplihed into helds on galactic scales. Considerable attention has been focused 
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on generation of seed fields during cosmic phase transitions, e.g., Ref. particularly at the 
electroweak transition [^], or at the QCD transition through formation of electromagnetic 
helds in the collision of shock fronts or formation of dipole charge layers on the surfaces of 
phase transition bubble walls whose motion generates an electromagnetic current and hence 
a magnetic held |^. Always present as a source of seed helds are random magnetic held 
huctuations on a scale of order a thermal wavelength. Once one has identihed the source 
of microscale seed helds, it is necessary to hnd how they are transformed into macroscopic 
scale magnetic helds, a process requiring production of current loops on large size scales, 
and thus the generation of a level of vorticity in the huid. 

In this paper we focus on the generation and amplihcation of helds at the electroweak 
phase transition, by means of the huid, and hence magnetic, turbulence which should be 
formed there. During a hrst-order electroweak phase transition bubbles of the new phase 
form 1^- [but see arguments in Ref. [1^ why the transition, even if hrst order transition, 
need not necessarily proceed dynamically via domain wall formation]. Then, as we argue in 
Sec. II, bubble collisions generate a level of turbulence and hence vorticity in the huid. As 
a bubble propagates through the electroweak plasma, it generates a precursor shock front 
which accelerates the huid outside the burning front of the bubble. As the shock fronts 
from diherent bubbles collide, they generate turbulent jets on a scale of order the size of 
a typical bubble. The turbulence in the huid amplihes whatever seed helds are present to 
hnite-amplitude large-scale size magnetic helds 0. These helds can be considerable, if 
there is time for equipartition between the kinetic energy associated with the turbulent how 
and the magnetic held energy. 

A possible mechanism of generation of magnetic seed helds is via the dipole electromag¬ 
netic charge layer on the surfaces of the bubbles that forms, as discussed in Sec. Ill, as a 
consequence of baryon asymmetry and the large mass of the top quark. The burning sur¬ 
faces of bubbles propagating within a region of turbulent how rotate in the moving huid. 
The rotation of the dipole charge layer thus sets up a current in the huid. We estimate the 
typical magnetic held generated by a single bubble, and then argue that these helds provide 
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at least one source of seed magnetic fields in the turbulent fluid which are then amplihed 
into full scale magnetic turbulence. The strength of the helds produced this way are consid¬ 
erably smaller in the absence of magnetic turbulence than the helds that would result from 
equipartition in the turbulent how. 

In order to estimate the generated magnetic held on scale sizes large compared to the 
electroweak scale through the cumulative ehect of the many bubbles in the system, we follow 
the classic analysis of Hogan |^, and determine, in Sec. IV, the magnetic held correlation 
function. 


C'(r) = (B(r)B(0)). (1) 

If magnetic helds dihuse up to a scale R, then the average magnetic held squared which 
remains after dihusion will be given in terms of the correlation function by 

(B^)^~{C(r))„, (2) 

where (■ ■ ■)^ denotes the spatial average over a volume of size ~ R?. From this result we 
estimate the typical energy in the magnetic held at any time at distance scales greater than 
the magnetic held dihusion length at that time. Finally we discuss the possible relevance of 
this work for cosmology. Throughout we use units in which c = % = 1. 


II. ELECTROWEAK PHASE TRANSITION BUBBLES AND TURBULENCE 


The electroweak phase transition appears, within acceptable parameters, to be weakly 
T^ . As the universe cools through the electroweak transition temperature. 


hrst order | 

Tc ~ 100 GeV, the plasma in the unbroken phase supercools; eventually, one assumes, small 
regions of the broken symmetry phase, with non-zero Higgs expectation value, nucleate, 
forming bubbles which expand and hll the system ||T^ - El- 

The typical size of a bubble after the phase transition is completed is in the range 


R 


bubble 




-1 

ew 


( 3 ) 
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where 


tr-i in 

hiew ~ ~ lUcm 


( 4 ) 


gl'X^ 

is the size of the event horizon at the electroweak scale, mpi is the Planck mass, ~ 10^ 
is the number of massless degrees of freedom in the matter, and the fractional size /b 
is ~ 10“^ — 10“^ These numbers are typical of parameters characteristic of the 

electroweak phase transition. 


The bubble wall surface appears to be stable against small fluctuations |]^ . The bubble 
wall velocity, while poorly known, should be in the range 


^^wall ~ 0.05 - 0.9. 


( 5 ) 


When Uwaii < 1/v^, the sound velocity in the symmetric phase of the electroweak plasma, 
the burning of the symmetric phase proceeds by deflagration in which the phase transition 
burning front expands out into the symmetric phase. Ahead of the burning front, a super¬ 
sonic shock moves into the symmetric phase, accelerating it outward, as shown schematically 
in Fig. la. On the other hand, when the bubble wall velocity is supersonic in the symmetric 
phase, the burning discontinuity generates in its wake a similarity rarefaction wave, with 
velocity and temperature prohles illustrated in Fig. lb. 

As two shock fronts associated with the nucleation bubbles collide, they generate tur¬ 
bulence, shown in Fig. 2. At the point where the shock fronts intersect, the velocity of 
the fluid from one bubble is in a different direction than that from the other bubble. The 
geometry is similar to the classic problem of the formation of a turbulent jet. The Reynolds 
number. Re, of the flow around two colliding bubbles is sufficiently large that we expect 
fully developed turbulence in a cone associated with the intersection of the two bubbles, 
whose opening angle depends on the angle of intersection of the two fluids ig. 

The Reynolds number for the collision of two bubbles is 


Re ~ 


^fluid -Rbubble 


A 


( 6 ) 
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where Ufluid is the typical fluid velocity, which we take to be ~ Uwaii ~ 10“^; the typical size 
of a bubble, -Rbubbie, is given by Eq. (H). The typical scattering length A of excitations in 
the plasma is of order 


A 


T luttwl 

where is the fine structure constant at the electroweak scale, and 


( 7 ) 

is the number 


of degrees of freedom that scatter by electroweak processes Therefore 


Re ~ 10 Inawl ~ 10^^ (8) 

is sufficiently huge, for any macroscopic size bubble, that the collision of the shocked matter 
should generate turbulent flow. 

A turbulent conducting fluid develops magnetic turbulence, resulting in magnetic fields 
on all scale sizes. The relevant time scale for the amplification of fields on length scale I is of 
order (V7?bubbie)ttrans, where ttrans ~ -Rbubbie/'Inwall E the duration of the phase transition. If 
the held growth is exponential, fields on scales /<i?bubbie can be amplified by many e-folds. 
When magnetic turbulence becomes fully developed, the kinetic energy of the turbulent flow 
is equipartitioned with that of the magnetic held energy, implying that the magnetic fields, 
-B(-Rbubbie), generated on the scale of the phase transition bubbles, the largest scale typical 
of the turbulent flow, are given by 


-B^(-Rbubbie) ~ e(Tc)nfl^id, 


(9) 


where e(Tc) ~ g^T^ is the energy density of the electroweak plasma. Since the velocity of the 
turbulent flow is ~ 10“^, a very large fraction of the energy of the fluid is in electromagnetic 
fields. Turbulent flow is therefore capable of amplifying any seed held by many orders of 
magnitude. 

Thermal fluctuations provide a omnipresent source of seed fields whose size can be es¬ 
timated by Eqs. (1) and (2), with the average in Eq. (1) the thermal average for a free 
photon gas. For scale R S> T~^ one has 
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( 10 ) 


i.e., for R ~ RhuhUe the value of -Bluet is smaller than the equipartition value by a factor of 
i9l^^/vlnid)fh^iTc/mpi)\ 

III. MAGNETIC FIELD GENERATED BY DIPOLE CHARGE LAYER IN THE 

BUBBLE WALL 

Let us now compare the magnitude of magnetic fields generated by turbulence with those 
generated from currents in the bubble walls in the absence of magnetic turbulence. These 
currents arise from the electric dipole layer that develops in the bubble wall, as a consequence 
of the baryon asymmetry of matter undergoing the transition combined with the fact that 
the mass of the top quark is comparable to the phase transition temperature. 

The net baryon number is non-zero in the broken symmetry phase. If baryogenesis is 
not associated with the electroweak phase transition, then the local baryon density in the 
neighborhood of the bnbble wall is characterized by a baryon number chemical potential 
Hi, ~ 10“®Tc. On the other hand, if the baryon asymmetry is generated at the electroweak 
phase transition, the asymmetry is driven by CP violating effects and /ij, can be mnch larger. 
We take as an acceptable range 10“® < iJ-b/R < 10“^. 

The expectation value of the Higgs field goes from zero to a finite value as one traverses 
the bubble wall into the broken symmetry phase. Thus the top quark has a non-zero mass 
inside the bnbble, and zero mass ontside. Top qnarks present ontside the bubbles face 
a potential barrier at the bnbble wall, which leads to a Boltzmann suppression of their 
nnmber inside the bubble and near the snrface. The slight excess of top over anti-top quarks 
contributes a net positive charge ontside the bnbbles, but the charge excess from the tops 
is suppressed inside. Overall charge nentrality is guaranteed by the presence of a small 
electrostatic potential, which is more negative inside the bnbble than outside, attracting 
light positively charged particles to the interior. 
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The electrostatic potential is determined by 


- = epem(a:), (11) 

where pem is the total local electric charge. Expanding the distributions of the various 
charged species to leading order in and the small non-zero charge, baryon and lepton 
number chemical potentials, and then for simplicity expanding in powers of (mtop/T)^, we 
hnd 


epem + -eml,p(x)lib, 


( 12 ) 


where mo is the Debye mass for massless leptons and quarks (including the top). (We 
neglect a small correction to the Debye screening term ~ The Debye screening 

length 1/mD is of order 1/eT. In the expected limit where the thickness, T^aii, of the bubble 
wall is S> 1/mD, we solve Eqs. (|TT]) and ([T^) by expanding in powers of V^/m^, and hnd 


/ \ f'h'fe ^2 2 / \ 

epem{x) -yV [x) 


ml 


(13) 


plus corrections involving higher order powers of l/(mDT„aii)^- The charge density on the 
surface, x ~ 0, is of order 


epem(O) ~ -VbT^e, 


(14) 


where the typical baryon asymmetry near the bubble wall is rji, Pfe/Tc, mf) and 

2 


^ _ j mtop 


( 16 ) 


\T)?L„all / 

Typically TcLwaii ~ 10-100, so that e ~ 10“^ — 10“^. For mtop{x) monotonically decreasing 
from inside to out across the bubble surface, we see that the bubble has a dipole charge 
layer of order epem(0)Twaii per unit area. 

In order for the dipole charge layer on the bubble wall to generate a magnetic held the 
bubble must have a net rotation. If the bubble is propagating in a turbulent region with 
vorticity, its surface acquires a net rotational velocity, Urot, of order the typical turbulent 
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velocity in the fluid, Ufiuid ~ 10 A rotating bubble thus has associated with it a magnetic 
moment of order 



(16) 


The magnetic held generated by the bubble on a scale size of order the bubble radius is 
therefore 



,2 Vb ^top 
" /b mompi ■ 


gl'^Vi^uT^ 


(17) 


We note that the magnetic held produced by this mechanism is smaller by a factor ~ 


rjfjTc/mpi compared with the held expected from equipartition in a magnetically turbulent 


environment, Eq. (^). At the electroweak scale the equipartition held is ~ 10^"^ gauss, and 
the held produced from the dipole layers in the bubble walls, ~ 10“^ gauss, is a factor 
10 smaller. 

Cheng and Olinto [0, by contrast. And a field produced by the dipole layer in the bubbles 
at the QCD transition, at T ~ 100 MeV, of order 10“^° — 10“^^ of the equipartition held, 
~ 10^® gauss. The reason for their much larger fraction is that they assume that the width 
of the dipole layer is controlled by the width of the baryon diffusion layer, ~ 10^ fm, rather 
than the microphysics length scale ~ 1 fm at the QCD scale. The diffusion layer for the 
baryon number is so much larger because the QCD transition occurs sufficiently slowly that 
the particles have time to undergo many scatterings during the transition. In the electroweak 
transition on the other hand, the bubble walls move relativistically, and little diffusion takes 
place. 


IV. LARGE SCALE MAGNETIC FIELDS 


We now ask how magnetic fields on scales larger than the electroweak bubble size are 


generated by superposition of the fields associated with the turbulent flow. A first estimate 
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can be made by assuming that the magnetic fields on scales large compared with the size of 
the turbulent velocity fields are generated by randomly oriented magnetic dipoles of typical 
size that of this velocity field. 

A magnetic dipole M generates a magnetic field far from the dipole of order 


g(f)~e , ^ „ (18) 

I r - Fd h 

where rd is the position of the dipole. Near the dipole the field saturates to a more or less 
constant value. To estimate the net magnetic fields produced by randomly oriented dipoles, 
we make a continuum approximation for the distribution of dipoles, and assume that the 
density ) of dipoles pointing in the ith direction is Gaussianly distributed with measure 


J [dz/*] exp 



(19) 


where k is a constant. The correlation function of the density of dipoles implied by the 
distribution (|I^) is 




( 20 ) 


From Eqs. 


and (pop we thus find the magnetic field-magnetic field correlation function 


(B(r) • B(0)) ~ 



( 21 ) 


This integral representation is dominated by the regions where either | rd 0 or | r —rd |—> 
0. The logarithmic divergence of the integral in these regions is cut off by the size of the 
typical dipole, fhH-J, so that for r > 

( 22 ) 

The average strength of measured by averaging on a size scale R is thusQ 


^This results differs from Hogan 


who on the basis of a random walk of the field lines finds 
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^ — In^ 

^ V /b . 


( 23 ) 


If we assume equipartition, and compare with Eq. (I), we see that k is given by 


e^K ~ 


(/^.^^ewh 




(24) 


and therefore 


2 \ 2 rr .4 f fh \ ,„2 / HewR\ 


R \ - Vflnid9*T, 


In" 


(25) 


V /b y 

The expansion of the universe causes the magnetic field to decrease as the square of scale 
factor of the universe. Therefore, in the absence of flux diffusion, the ratio of {B^) to the 
energy in photons. 


1?B = 


P-t 


(26) 


remains independent of time, when measured on a comoving scale size. From (^^ we have 
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(27) 


where Aew is the size of the universe at the electroweak phase transition times the scale 
factor, Tc/Ty, where is the temperature of the microwave background; at the present 
epoch Aew ~ 10^ AU. 

Magnetic fields from the electroweak transition can survive only on scales on which 
magnetic diffusion has not had time to wash out the field correlations. The flux diffusion 
equation, which in the local fluid rest frame can be written as. 


+ AiradtB = 0 , 


(28) 


where a is the electrical conductivity, implies that a magnetic held of length scale 1/k dies 
off on a time scale 


T ~ cr/k'^; 


(29) 
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long wavelength magnetic fields set up at the electroweak phase transition die away very 
slowly. 

The characteristic diffusion distance at time r is therefore of order 



(30) 


at T magnetic fields on length scales larger than are not yet significantly affected by 


diffusion. An upper limit on the diffusion length is its present value, about 3 AU 0] which 


is not much smaller than the size of the horizon at the electroweak phase transition. 

In the present epoch, rjB ~ (10“® — 10“®) x (lO^AU/i?)^ In^(i?/1AU) for ft, = 10“^ — 
10“^. Thus the equipartition magnetic field on a scale of order of diffusion length in the 
extragalactic medium is now about ~ lOAU) ~ 10“^ — 10“® gauss. These fields might 
have acted in an earlier epoch as seed fields, subsequently ampified by a galactic dynamo 
mechanism to produce the galactic magnetic fields. On a galactic size scale, at present seven 
orders of magnitude larger than the present size of the electroweak horizon, one finds that 
the equipartition field is ~ 10®AU) ~ 10“^^ — 10“^® gauss. 
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FIGURES 

FIG. 1. Bubble structure in (1+1) dimensions for a first-order electroweak phase transition. 
The energy density e and the fluid rapidity 0 = ln[(l + Ufiuid)/(1 — ffluid)] are plotted versus 
the space-time rapidity y = ln[(t + x)/{t — x)]. (a) Deflagration bubble; 0def and 0sh denote 
the rapidities of the deflagration front and the shock front, respectively, (b) Detonation bubble; 
here 0s is the rapidity corresponding to the sound velocity and 0det denotes the rapidity of the 
detonation front. 

FIG. 2. Fluid velocities during the collision of two bubbles. The dashed line represents the 
shock (detonation) front for a deflagration (detonation) bubble. 
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